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Abstract 

The Langley 8-Foot High- Temperature Tunnel 
(8-ft HTT) is used to test components of hyper- 
sonic vehicles for aerothermal loads definition and 
structural component verification. The test medium 
of the 8-ft HTT is obtained by burning a mixture 
of methane-air under high pressure; the combus- 
tion products are expanded through an axisymmet- 
ric conical-contoured nozzle to simulate atmospheric 
flight at Mach 7. This facility has been modi- 
fied to raise the oxygen content of the test medium 
to match that of air and to include Mach 4 and 
Mach 5 capabilities. These modifications will facil- 
itate the testing of hypersonic air-breathing propul- 
sion systems for a wide range of flight conditions. 
A computational method to predict the thermo- 
dynamic, transport, and flow properties of the 
equilibrium chemically reacting oxygen-enriched 
methane-air combustion products has been imple- 
mented in a computer code. This code calculates the 
fuel, air, and oxygen mass flow rates and test sec- 
tion flow properties for Mach 7, Mach 5, and Mach 4 
nozzle configurations for given combustor and mixer 
conditions. Salient features of the 8-ft HTT are de- 
scribed, and some of the predicted tunnel operational 
characteristics are presented in the carpet plots to as- 
sist users in preparing test plans. 

Introduction 

The Langley 8-Foot High-Temperature Tunnel 
(8-ft HTT), which became operational in the 1960's, 
primarily has been used to define thermal and struc- 
tural loads on large models at hypersonic speeds. 
Although many other hypersonic facilities were aban- 
doned in the early 1970’s, this tunnel is still opera- 
tional and qualifies as a unique national facility. A 
resurgent interest exists in the development of air- 
breathing hypersonic vehicles that can fly at hyper- 
sonic speeds within the atmosphere. This interest 
is evident from the current development program of 
the National Aero-Space Plane (NASP), which is en- 
visioned to be a hypersonic aircraft that can fly di- 
rectly into orbit from a conventional runway (refs. 1 
and 2). 

To develop and improve advanced propulsion sys- 
tems required for vehicles such as NASP, test fa- 
cilities must simulate conditions that are encoun- 
tered by the vehicle during its mission. Although 
the high-temperature test medium of the methane- 
air combustion products produced in the 8-ft HTT 
does not have the oxygen content of air, it has been 
useful for aerothermal loads definition and structural 
component verification. However, engine testing re- 
quires oxygen enrichment because most of the oxy- 


gen available in the air is used in the combustion of 
methane and the test medium cannot support fur- 
ther combustion in the engine. The need for large- 
engine test facilities to develop air-breathing engines 
for hypersonic flight prompted the modification of 
the 8-ft HTT. The modified facility can be operated 
in the oxygen-enriched combustion mode to produce 
21 percent oxygen by volume in the test section or 
in the no-oxygen-enriched mode. The addition of 
Mach 4 and Mach 5 nozzle configurations will also 
complement the existing Mach 7 capability. 

Figure 1 illustrates the operational envelopes, in 
terms of Mach number and pressure altitude, for the 
8-ft HTT and two other facilities (the Aero Propul- 
sion Test Unit (APTU) and the Aero Propulsion 
System Test Facility (ASTF) of the U.S. Air Force 
Arnold Engineering Development Center (AEDC)). 
Also superimposed around the facility operational 
envelope is a typical air-breathing engine operational 
envelope. The upper altitude limit is imposed by the 
ability of the engine to sustain combustion, and the 
lower altitude limit is imposed by the ability of the 
structure to survive the aerothermal loads. The in- 
clusion of Mach 4 and Mach 5 capabilities enhances 
the testing envelope of the 8-ft HTT in the turbo- 
jet and ramjet operating ranges. However, the lower 
altitude limit of the facility with oxygen enrichment 
is restricted because the LOX (liquid oxygen) tank 
pressure limit is 2300 psia, which limits the combus- 
tor pressure to 2000 psia. 

This report describes a computational method to 
predict test section flow properties at the nominal 
Mach numbers of 7, 5, and 4 for the 8-ft HTT. 
The flow is assumed to be in chemical equilibrium. 
Salient features of the 8-ft HTT are discussed to make 
the code requirements clear. The tunnel operational 
characteristics are presented for both the methane- 
air and methane-air-oxygen modes. 


Symbols 


A 

nozzle throat cross-sectional area 

a 2 

a. 

acoustic velocity, ft /sec 

c p 

molar heat capacity of species, 
Btu/molc-°R 

Cp 

specific heat capacity of gaseous 
mixture, Btu/lbm-°R 

FSA 

flow survey apparatus 

9 

acceleration due to gravity, 
ft /sec 2 


H 

h 

J 

K 


K, 


p.j 


LOX 

M 

m 

MW 

ATpr 

N Ke 

P 

PTC 

PTM 

Q 

9 

9s 

R 

Rn 

Rs P 

S 

s 

T 

TTC 

TTM 

V 


sensible enthalpy of species, 
Btu/mole 

enthalpy of gaseous mixture, 
Btu/lbm 

mechanical equivalent of heat, 
778.26 ft-lbf/Btu 

thermal conductivity, Btu/ 
ft-sec-°R 

equilibrium constant in terms 
of partial pressures (eqs. (A7) 
to (A12)) 

Boltzmann constant, Btu/ 
molecule-°R 

liquid oxygen 

Mach number 

mass flow rate, lbm/sec 

average molecular weight of 
gaseous mixture, Ibm/mole, 1/a 

Prandtl number 

unit Reynolds number, ft -1 

pressure, psia 

combustor pressure, psia 

mixer pressure, psia 

heat of formation, Btu/mole 

dynamic pressure, psi 

stagnation point heat-transfer 
rate, Btu/ft 2 -sec 

universal gas constant, 

1.9858 Btu/mole-°R 

radius of spherically blunt body, 
1 ft 

specific gas constant for mixture, 
Rj MW, Btu/lbm-°R 

entropy of species, Btu/mole-°R 

entropy of gaseous mixture, 
Btu/lbm-°R 

temperature, °R 

combustor temperature, °R 

mixer temperature, °R 

velocity, ft /sec 


x proportion of oxygen in gaseous 

mixture by volume 

y moles of oxygen per mole of fuel 

z moles of air per mole of fuel 

7 ratio of specific heats 

e convergence criteria 

A Lennard-Jones force constant 

fi viscosity of gaseous mixture, 

lbm/ft-sec 

( collision diameter of molecule, A 

p mass density of mixture, lbm/ft 3 

a mole number of mixture, 

mole/lbm of mixture, E 

moles of elements C, H, O, and N 
per lbm of mixture, mole/lbm 


Q 

Subscripts: 


e 

/ 


771 

OX 

r 

s 

t 

Uj d 
1,2 

Superscript: 


collision integral 
air 

combustor air 

effective 

fuel 

chemical species index; pre- 
combustion species 

chemical reaction index; post- 
combustion species 

mixer air 
molecular oxygen 
reference 

sum of fuel, total air, and oxygen 
total 

upstream and downstream 
sections of transpiration-cooled 
nozzle (fig. 8) 

locations shown in figure 8 
computed value 


Test Facility 

General Description 

A schematic of the 8- ft HTT is given in figure 2. 
The tunnel is a large, hypersonic blowdown facility 
that simulates true- temperature flight at altitudes 
and Mach numbers shown in figure 1. The facility 
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obtains its high-energy test medium by burning a 
mixture of methane and air under high pressure 
in the combustor. The combustion products are 
expanded to the test stream Mach number through 
an axisymmetric conical-contoured nozzle that has 
an exit diameter of 8 ft. The free- jet flow passes 
through the test chamber and enters a straight-tube 
supersonic diffuser. The flow then is pumped by a 
single-stage annular air ejector into a mixing tube 
and exhausted to the atmosphere through a subsonic 
diffuser. 

A view of the tunnel test section from the nozzle 
exit during the entry of a model into the test stream 
is shown in figure 3. The model is stored in a pod 
beneath the test stream to protect it from tunnel 
start-up loads. Once the flow conditions are estab- 
lished, the model is inserted into the stream using 
a hydraulically actuated elevator. Model insertion 
time from the edge of the test core to the tunnel 
centerline is approximately 1 sec. An array of quartz 
lamps, located in the pod, is used to radiantly heat 
the models, if required, before and after exposure of 
the model to the test stream. A flow survey appa- 
ratus (FSA), which can obtain measurements of the 
test section flow properties and can survey the flow 
at various axial locations along the test section, is 
shown in its stowed position. This FSA is instru- 
mented with a total of 37 probes, including 13 pitot, 
11 static pressure, and 13 total- temperature probes. 

The combustor operates at total temperatures 
ranging from 2400°R to 3600°R. For the Mach 7 
nozzle configuration with no oxygen enrichment, the 
free-stream dynamic pressure ranges from 1.74 psi to 
12.5 psi, and the free-stream unit Reynolds number 
ranges from 0.3 x 10 6 per ft to 2.2 x 10 6 per ft. The 
maximum run time of the tunnel is approximately 
120 sec, and it is primarily restricted by the air 
storage capacity of a local bottle field that is the 
common source for the combustor and the ejector. 

Figure 1 shows that the test conditions for the 
Mach 7 nozzle configuration are at the lower end 
of the scramjet range. To obtain test conditions 
for lower Mach numbers and pressure altitudes, the 
nozzle must be reconfigured. A dual-throat-mixer 
concept similar to that used to convert the AEDC 
Von Karman Gas Dynamics Facility from a Mach 10 
tunnel to a Mach 4 true-temperature tunnel (ref. 3) 
was employed to reduce the cost by utilizing the max- 
imum portion of the existing facility. (Note that the 
Mach 10 tunnel is not a true-temperature tunnel.) 
However, in the present facility, gas dynamics rather 
than mechanical means (ref. 3) are used to mix the 
hot combustion products with air. A similar concept 


has been used to provide increased capabilities for the 
Langley Arc-Heated Scramjet Test Facility (ref. 4). 

The actual Mach number for the Mach 7 nozzle 
configuration varies from 5.8 to 7.3 (fig. 1). This 
Mach number variation is caused by the water con- 
densation in the nozzle (ref. 5). Condensation is 
found to be maximum at low temperatures and high 
pressures and minimum at high temperatures and 
low pressures. In the future, the combustor tem- 
perature upper limit of 3600°R may be increased 
to 4000°R with an improved thermal liner for the 
combustor. 

Oxygen Enrichment 

The 8-ft HTT has been modified to include oxy- 
gen enrichment for facilitating scramjet engine test- 
ing. A schematic of the combustor, which shows the 
injection location and path of the LOX along with the 
fuel and air supply, is given in figure 4. High-pressure 
air from the 6000 psi bottle field is introduced at 
the upstream end of the combustor. The combus- 
tor consists of a carbon steel outer pressure vessel 
that is protected by a 316 Stainless-Steel outer liner 
and a 201 Nickel inner liner. The supply air passes 
downstream through the annulus between the pres- 
sure vessel and stainless-steel liner, and then it turns 
180° upstream in the annulus between the outer and 
inner liners. The LOX enters the annulus between 
the outer and inner liners near the end of the LOX 
injection mixing ramp and mixes with the upstream 
air before it exits the annulus. The flow once again 
turns 180° and passes across the methane injectors 
where the oxygen-enriched air mixes with methane 
and burns. The combustion gas then flows through 
the nozzle throat and expands into the test section. 

Because the LOX run tank pressure is limited 
to 2290 psia, the combustor pressure is limited to 
2000 psia. Therefore, the facility that operates with 
oxygen enrichment is limited to a maximum dynamic 
pressure of approximately 1800 psf or minimum pres- 
sure altitudes of about 60 000 ft at M = 4, 78 000 ft 
at M = 5, and 90 000 ft at M = 7. 

Air- Transpiration-Cooled Nozzle Throat 

Prior to the facility modification, the nozzle 
throat was air-film cooled to prevent it from over- 
heating because of the high- temperature combustion 
products. The water-cooled nozzle approach section, 
the throat, and a portion of the nozzle expansion 
section have been replaced with an air-transpiration- 
cooled section (fig. 5). This section is approximately 
7.2 ft long and has a maximum internal diameter of 
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36 in. and a throat diameter of 5.6 in. A platelet 
concept has been used to fabricate this section in 
which the air passages are photoetched on thin sheets 
of varying thicknesses and bonded together. This 
air-transpiration-cooling approach is expected to re- 
duce the coolant air mass flow rate by 50 percent of 
the amount that was previously used just to air-film 
cool the throat. This reduction in the coolant flow 
should produce a larger, uniform-temperature test 
section core. The coolant flows at the upstream and 
downstream locations of the nozzle throat are con- 
trolled such that they are proportional to the com- 
bustor pressure. However, the optimum values of the 
coolant flows will be determined experimentally so 
that the boundary layer does not separate from the 
nozzle wall but cools the nozzle effectively. Refer- 
ence 6 gives more information concerning the design 
details of the air- transpiration-cooled section. 

Alternate Mach Number 

Schematics of the nozzle modifications to include 
Mach 4 and Mach 5 capabilities are shown in fig- 
ure 6. As explained in the section entitled ‘'General 
Description,” a dual-throat-mixer concept has been 
selected to preclude modifications to the combustor 
and maintain the 8- ft nozzle exit diameter for Mach 4 
or Mach 5 operation. For these Mach numbers, the 
removable section of the Mach 7 nozzle is replaced 
with a mixer and the Mach 4 or Mach 5 nozzle in- 
sert that joins the 18.4-ft-long fixed nozzle section 
as shown in figure 6. Air at ambient temperature 
is injected at various locations of the mixer to de- 
crease the temperature and momentum of the hot 
gas and increase the mass flux in the mixer. The gas 
then will be expanded through the Mach 4 or Mach 5 
nozzle to simulate the true temperature of flight at 
realistic altitudes. The gas temperature in the mixer 
is 1640°R and 2350°R for the Mach 4 and Mach 5 
configurations, respectively. 

Flow Computational Method 

General Description of Code 

The computational sequence of the computer code 
to calculate the thermodynamic, transport, and flow 
properties for the 8- ft HTT is described in the flow- 
chart in figure 7. This code, called HTT, calculates 
the test stream flow properties and the mass flow 
rates of fuel, air, and oxygen at various locations of 
the 8- ft HTT for the Mach 7, Mach 5, and Mach 4 
nozzle configurations (fig. 8). These mass flow rates 
will be monitored and controlled during the tunnel 
operation. 

The combustion products are assumed to be in 
chemical equilibrium and to be composed of 4 ele- 


ments (C, H, O, and N) and 10 reacting species (H20, 
C0 2 , CO, 0 2 , H 2 , N 2 , H, O, OH, and NO) which 
obey the perfect gas law. First, the proportions of 
fuel, air, and oxygen are determined using an iter- 
ative procedure such that the combustion products 
attain the specified percentage of oxygen by volume 
and temperature at a given pressure. Note that the 
methane-air mode is calculated by setting the oxy- 
gen volume to the amount contained in air. Calcu- 
lating the temperature of the combustion products 
at a specified pressure requires the mole numbers 
<r c ? cqp (Jo? and a ^ of the gas supply. These mole 
numbers are used to determine the chemical com- 
position by solving a set of simultaneous equations 
which represent the chemical reaction and elemen- 
tal mass balance. The thermodynamic properties of 
the mixture are then computed using thermodynamic 
equations and chemical tables. The temperature of 
the combustion products is calculated assuming that 
the enthalpy of post combust ion and precombustion 
products is constant for a given combustor condition. 
The proportion of the gas supply is adjusted until the 
calculated temperature converges with the specified 
temperature. The actual quantities of fuel, air, and 
oxygen are then calculated. 

Next, the test stream flow properties are com- 
puted using either the FSA data or isentropic ex- 
pansion procediire. In the first method, the values 
of test stream static pressure p\ and pitot, pressures 
pt 2 are obtained from the FSA data. These values 
are used in the governing equations to obtain all the 
other flow properties. The flow properties for the 
Mach 7 nozzle configuration are computed using the 
FSA data, which are available from measurements 
obtained prior to the tunnel modification. Calcu- 
lations for the Mach 4 and Mach 5 nozzle configura- 
tions are done assuming isentropic expansion because 
the FSA data are not presently available. For this 
method, p\ and p^ 2 are calculated by first isentrop- 
ically expanding the gas to the desired test stream 
Mach number and then solving the governing equa- 
tions. Once the values for p\ and p^ 2 are obtained, 
the calculation procedure for all the other test stream 
flow properties is the same for both these methods. 
Finally, the transport properties of the mixture are 
calculated. 

Combustor Conditions 

Gas mixture . The calculation method for the 
proportions of fuel, air, and oxygen enrichment to 
produce the combustion products at specified com- 
bustor pressure PTC and combustor temperature 
TTC with the desired volumetric proportion of oxy- 
gen x is illustrated by the initial portion of the 
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computer flowchart in figure 7. The chemical equa- 
tion for methane-air with oxygen enriched combus- 
tion is 

CH 4 + z(0.21O‘2 + 0.79N 2 ) + y0 2 

= C0 2 + 2H 2 0 + {y + 0.21* - 2)0 2 + (0.79z)N 2 (1) 

This equation represents the chemical reaction for 
1 mole of methane fuel with z moles of air and y moles 
of oxygen enrichment. The air is assumed to have 
a nitrogen to oxygen volumetric ratio of 79 to 21. 
The numerical coefficient of O 2 , y 4- 0.21z — 2, in 
the right-hand side of the equation corresponds to 
the amount of oxygen present in the combustion 
products. 

The proportions of oxygen and air in terms of 
methane fuel can be used to define the elemental 
composition of <7c, an d respectively. For 

CH 4 

^h/^C = ^ (2) 

The mole ratios of oxygen to carbon and nitrogen to 
carbon from the left-hand side of equation (1) are 

*o/*C = 2(y + 0.21s) (3) 

and 

^n/^C = 2(0. 79z) (4) 

The four elemental constants of <jq , <jq, <th> 
and cqv per unit mass of mixture are related by the 
following elemental mass balance equation: 

12ctq *T <th + 16ao T 14 <tn = 1 (5) 

where the numerical coefficients of oq, crjp 
and <tn are the atomic weights of C, H, O, and N, 
respectively. An expression for ctq is obtained by 
combining equations (2) to (4) with equation (5) as 

<jq = 1/(16 T 32y -|- 28.84^) (6) 

The proportion of each clement supplied to the 
combustor is a function of moles of oxygen and air 
per mole of methane fuel. This proportion can be 
resolved by specifying the volumetric ratio x c of the 
oxygen to the total of the products in the combustor 
defined by the following equation: 

x c = {y + 0.21z c - 2 )/(y + z c + 1) (7) 

The numerator of this ratio is the coefficient of 
oxygen in the combustor which is contained in the 
third term of the right-hand side of the general 
equation (1) with z c instead of z. The denominator 


is the total moles of the reacting combustor gas in 
the equation. In general, x and z arc defined to 
describe the gas mixture that is produced in the 
test stream. The subscripted parameters x c and 
z c (fig- 7) correspond to the combustor gas that is 
directly expanded to the test stream for the Mach 7 
case or the mixer for the Mach 4 and Mach 5 cases. 

As indicated by the flowchart in figure 7, the 
gas composition supplied to the combustor is deter- 
mined by specifying the required proportion of oxy- 
gen in the test section x (which is equal to x c for the 
Mach 7 case) and by assuming an initial value for z c 
to compute y from equation (7). The combustor el- 
emental mole numbers are computed using the com- 
bustor air moles z c instead of z in the general equa- 
tions (2) to (4) and (6). The chemical composition 
of the gas is then computed as described in appen- 
dix A; 10 species and 6 chemical reactions (eqs. (Al) 
to (A6)) are considered, and the species are assumed 
to be in chemical equilibrium. The chemical compo- 
sition is computed for a specified combustor pressure 
and temperature by simultaneously solving the set 
of equations (A7) to (A 16) (see ref. 7 for details). 
The temperature of the combustion products is then 
computed as explained in appendix B. Iterations are 
made on z c until the computed temperature of the 
combustion products converges on the specified com- 
bustor temperature (fig. 7). This iterative procedure 
wall give the correct values of z c and y and the corre- 
sponding chemical and thermodynamic properties of 
the combustor gas. For the methane-air case, no oxy- 
gen enrichment exists, and the computation is done 
with y set to zero. 

Mass flow rates . The transpiration-cooled air 
is injected upstream and downstream of the Mach 7 
nozzle throat, respectively. The coolant mass flow 
rates are controlled at the designed level (which 
is proportional to PTC) and are defined by the 
respective equations: 

m u = 0.021(PTC) (8) 

m d = 0.017(PTC) (9) 

These correlations may be modified to include TTC 
based upon the future optimization analysis of the 
nozzle cooling performance. 

The mass flow rates of combustor air m c and oxy- 
gen m ox per unit mass flow rate of fuel are deter- 
mined from the computed values of z c and respec- 

tively, by using the following equations: 

mc/mj = z c (MW a /MW f ) (10) 
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I if min 


m ox /m f = y(MW ox/MW /) (11) The resulting mixture pressure is computed, assum- 

ing the second throat is choked, as 

However, to compute the actual mass flow rates of 
each gas, another equation that relates these mass 
flow rates is developed assuming that the Mach 7 
nozzle is choked. 


\m s 


\j (TTM) 


(0.5324) (14) 


First, a nozzle throat annular area that corre- 
sponds to a choked m u (eq. (8)) is computed assum- 
ing a total pressure of PTC and a total temperature 
of 1360°R, which is the nozzle design inner surface 
temperature. Next, an effective throat area, which 
is the Mach 7 nozzle geometrical throat area minus 
the annular area that corresponds to the coolant flow 
m u , is calculated. Finally, the choked mass flow rate 
through the effective area m e is calculated using PTC 
and TTC of the combustor gases. This mass flow rate 
m e is the sum of fuel, air, and oxygen which are sup- 
plied to the combustor, as indicated by the following 
expression: 

m e = mj + m c + m ox (12) 


The value of Zt is adjusted, and the gas mixture com- 
position and properties are computed until conver- 
gence on TTM is obtained utilizing the same proce- 
dure used for the combustor. The value of y for the 
mixer calculations corresponds to the LOX supplied 
to the combustor. 

For normal operations in which the desired oxy- 
gen content of the test stream x is 0.21, the combus- 
tor gas produced at x c = 0.21 is mixed with air (with 
the same oxygen content) so that x f = x. However, 
if a different oxygen content is desired, an additional 
iteration is required to adjust x C) and the whole com- 
putational procedure is repeated using zt instead of 
z c in equation (7) until the desired oxygen content is 
obtained (fig. 7). 


The mass flow rates of mf, m c , and m ox are calcu- 
lated by solving equations (10) to (12). 

Mixer Conditions 

For Mach 4 and Mach 5 nozzle configurations, 
the mixer supplies additional air to achieve flight- 
temperature simulation and higher mass flow rates 
for low- altitude simulation. The airflow rate to 
the mixer m m and pressure PTM are calculated to 
provide a specified mixer temperature TTM. The 
gas temperatures in the mixer for the Mach 4 and 
Mach 5 cases were assumed to be the design values 
of 1640°R and 2350°R (for flight-temperature sim- 
ulation), respectively. These temperatures produce 
a test stream static temperature of approximately 
400°R at the nozzle exit. However, the code can be 
run for a range of mixer temperatures. 

The computation procedure for the mixer flow 
quantities (fig. 7) is similar to that explained in the 
section entitled “Mass Flow Rates” for computing 
the combustor flow quantities. For the combustor 
analysis, z c moles are represented by m c , but for the 
mixer, z% moles of air are represented by the total 
airflow m^, which is the sum of combustor air ra c , 
the transpiration-coolant air m u and and the 
mixer air m m . The total mass flow rate through the 
second throat is the sum of the fuel and the oxygen 
supplied in the combustor and the total airflow rate, 
which is 

(13) 


Test Stream Properties 

As indicated by the latter portion of the computer 
flowchart in figure 7, the test stream properties are 
computed from data correlations measured using the 
tunnel FSA or assuming isentropic expansion of the 
combustion products to the test stream Mach num- 
ber. When using FSA data, p\ and pt } 2i respectively, 
are linearly correlated to TTC and PTC. In the latter 
case, pi and pt£ are calculated by isentropically ex- 
panding the gas to the desired Mach number and by 
solving the continuity, momentum, energy, and state 
equations. The other test stream properties, on both 
sides of the normal shock, are computed by assuming 
that the enthalpies ht,\ and ht , 2 are equal and by solv- 
ing the governing equations. The chemical compo- 
sition, thermodynamic properties, and other stream 
properties are calculated by methods described in ap- 
pendices A, B, and C, respectively. 

Flow survey apparatus data or isentropic 
expansion . The test stream properties, which are 
computed using FSA measurements of p\ and ptp 
(fig. 8), include the losses in the tunnel. The mea- 
sured pt , 2 and pi, at a location 84 in. downstream of 
the Mach 7 nozzle exit, are correlated to TTC and 
PTC as 

Pt, 2 /PTC = (1.2856 x 10" 6 )TTC + 0.00277 (15) 

and 

P\/Pt,2 = (-4.00 x 10~ 6 )TTC + 0.03080 (16) 


6 


m s =mj + m ox + rn t 


These correlations for the Mach 7 nozzle configura- 
tion, which are from FSA data that were obtained 
before the tunnel modification, will be corrected, if 
necessary, once FSA data from the modified tunnel 
are available. 

In the case of the isentropic expansion procedure, 
p\ and ptp are calculated by isentropically expanding 
the combustion products to the required test stream 
Mach number and by solving the mass, momentum, 
energy equations, and state equations. Because of 
the unavailability of FSA data for Mach 5 and Mach 4 
nozzle configurations, the isentropic expansion pro- 
cedure is used to compute the values of p\ and pt } 2- 
However, when FSA data become available, the cor- 
relation for p\ and pt } 2 can be expressed in terms of 
TTM and PTM, similar to equations (15) and (16), 
and can be incorporated into the code. From these 
computed pi and pt ? 2 values, the calculation method 
for all the flow properties upstream and downstream 
of the shock is described in the following sections. 

Computation of flow properties . The value 
of pt$ is obtained from equation (15) for given com- 
bustor conditions. The temperature of the combus- 
tion products Tt } 2 at 2 is calculated as explained in 

appendix B. The gas density at the stagnation con- 
ditions of pt$ and T t) 2 is computed from the state 
equation. The gas is isentropically expanded from 
p t 2 and T* 2 to P2 at an initially assumed value of T 2 
(fig. 7), and all the gas properties at P 2 and T 2 are 
determined as explained subsequently. 

First, a new static pressure p f 2 is calculated as 
P2=P2- dp 

where the incremental pressure dp is calculated from 
dp = -P 2 (s 2 ~ s t) 2 )/(*kr) 

The term p 2 is substituted for P 2 , and the iteration 
is continued until dp satisfies a convergence criteria. 
The gas composition and the thermodynamic prop- 
erties are calculated at P 2 and T 2 , and the flow prop- 
erties p 2 and V 2 are calculated from 

92 = P2/R&T2 (17) 

and 

V 2 = [2J(h t ,2 ~ h 2 )} 1/2 (18) 

Next, Vi and p\ y which are obtained by solving 
the mass and momentum equations across a normal 
shock, are given by 

V\ = \{p2 -Pl)/P2^2] + V 2 (19) 


and 

Pi = P2(V 2 /Vi) (20) 

The value of T\ is calculated from the state equa- 
tion using the values pi and pi obtained from equa- 
tions (16) and (20), respectively, and an initial as- 
sumption for cr. The chemical composition and a 
new value of cr are computed at pi and 7\, and T\ 
is updated until it converges on o. The enthalpy of 
gaseous mixture hi at pi and T\ is determined, and 
a differential enthalpy is calculated as 

dh — ht 2 ~ + 0.5V] 2 / 

The initially assumed T 2 in the outermost loop of this 
subroutine is adjusted, and the entire computation is 
continued until the value of dh is less than a conver- 
gence criteria e (fig. 7). The chemical composition, 
thermodynamic properties, and transport properties 
for the converged solution are calculated as described 
in appendices A, B, and C, respectively. 

The Mach number, unit Reynolds number, dy- 
namic pressure, and stagnation heat-transfer rate are 
calculated. The Mach number is defined as 

M = V/a (21) 

where a — y/g(dp/dp). In this equation, dp is calcu- 
lated as the difference of static pressures correspond- 
ing to two temperatures of 7* — dT and Tt + dT at pt. 
The value for dT is arbitrarily chosen to be 1°R. The 
difference in the densities dp, corresponding to dp for 
this change in Tt , is determined from the calculated 
chemical composition and the state equation. 

The unit Reynolds number is calculated from 

Afae = pV/p (22) 

The dynamic pressure is calculated from 

q = (0.5pV r2 )/144 (23) 

The stagnation heat-transfer rate q s for a spherically 
blunted body of a unit radius of 1 ft is calculated 
using the Fay and Riddell equation (ref. 8). The 
quantities of Cp, h, s, p, and p, which are required 
for the calculation of q s , are computed for a wall 
temperature of 540° R. 

Results and Discussion 

Mach 7 Configuration 

The carpet plots of the flow properties for the 
21 -percent oxygen enrichment and no-enrichment 
(no-LOX) cases are shown in figure 9 for the Mach 7 
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nozzle configuration. The combustor operating 
ranges for PTC and TTC are 600 psia to 3500 psia 
and 2500°R to 4000°R, respectively. However, for 
the oxygen enrichment case (shown by the dashed 
lines), PTC is limited to 2000 psia because of the 
oxygen run tank pressure limit. In the cross plots 
of pi versus T\ (fig. 9(a)), T\ is slightly higher for 
the oxygen enrichment case compared with the no- 
LOX case. This difference increases with increasing 
TTC. The variation of T\ with combustor pressure 
is insignificant for both the LOX and no-LOX cases. 
Note from equation (16) that p\ is assumed to be 
the same for both the LOX and no-LOX cases, and 
it will be updated, if necessary, once the tunnel sur- 
vey measurements are available. The difference in the 
properties of Vi, pi, M\, Np e , q s , and q \ , shown in 
figures 9(b) to 9(d), for the LOX and no-LOX cases 
is very small. The Mach number variation of approx- 
imately 6.1 to 7.2 in figure 9(c) is because of water 
condensation. This variation is accounted for by the 
FSA measurements in the test section. 

The test stream gas composition, in terms of mole 
fraction cr^/cr, is given in table I for both the LOX 
and no-LOX cases for the 21 -percent oxygen enrich- 
ment. Note that only the major species of H2O, 
C0 2 , 0 2 , and N 2 are listed in this table. The cal- 
culated free-stream static temperature is less than 
the chosen reference temperature of T r = 536.67°R, 
and the minor species are not considered for tem- 
peratures below T r (appendix A). The gas composi- 
tion is listed for the combustor design temperature 
of TTC = 3560°R. The composition is affected by 
temperature, but it docs not show significant change 
within the operating range of the combustor pressure 
of 600 psia to 3500 psia. The calculated value of the 
test stream gas constant R s p at this design temper- 
ature, for both the LOX and no-LOX cases for the 
21-percent oxygen enrichment, is 0.070 Btu/lbm-°R, 
and the corresponding values of the ratio of specific 
heats 7 are 1.38, 

Mach 5 and Mach 4 Configurations 

The carpet plots for the Mach 5 and Mach 4 noz- 
zle configurations for a combustor pressure range of 
600 psia to 3500 psia and combustor temperature 
of 3560°R are shown in figures 10 and 11, respec- 
tively. The mixer design temperatures are 2350° R 
and 1640°R for the Mach 5 and Mach 4 nozzle con- 
figurations, respectively. The plots are shown for the 
mixer design temperatures and for a few more mixer 
temperatures. A very small variation exists in the 
flow properties between the LOX and no-LOX cases. 
The variation observed in the plot p\ versus T\ for the 
Mach 7 case (fig. 9(a)) is not seen in similar plots for 


the Mach 5 and Mach 4 cases (figs. 10(a) and 11(a)), 
which is mainly because of the dominant effect of air 
added in the mixer. 

Figures 12 and 13 show the test stream properties 
in terms of PTM and TTM for the Mach 5 and 
Mach 4 cases, respectively. The plots in these figures 
can be used to estimate the test conditions from the 
mixer conditions alone. Correlations of combustor 
and mixer conditions for the Mach 5 and Mach 4 
cases are shown in figures 14 and 15, respectively. 
Note that no significant variation exists in these 
plots between the LOX and no-LOX cases. Some 
of the tunnel operational range and restrictions can 
be understood from these plots. For example, at 
TTM — 2350°R (the temperature required for true- 
temperature flight simulation) and PTM = 300 psia, 
which is shown by the bold line in figure 14(b), 
the combustor can be operated at different total 
temperatures. However, for this mixer condition and 
for TTC = 3000° R, the combustor pressure is above 
2000 psia (fig. 14(a)). The combustor cannot be 
operated at this temperature with the LOX because 
the LOX run tank pressure upper limit is 2000 psia. 
Figure 14(b) also shows that for the Mach 5 case at 
TTM = 2350° R, the combustor cannot be operated 
at TTC = 2500° R because the transpiration coolant 
flows (m u and m j in fig. 8) cool the combustor gas 
slightly below TTM = 2350°R. As a result, the code 
computes negative values of the mixer airflow m m , 
as shown in the plot in figure 14(b). However, the 
mixer will require air addition for cooling. 

Typical test stream gas composition for the 
Mach 5 and Mach 4 cases is given in table I at the 
mixer design temperatures. The gas composition is 
not affected by PTC and TTC, but it changes with 
TTM, Similar to the Mach 7 case, the minor species 
are set to zero for the Mach 5 and Mach 4 cases be- 
cause the test stream static temperature is below the 
reference temperature of T r = 536.67°R. The values 
of R s p are 0.070 Btu/lbm-°R and 0.069 Btu/lbm-°R, 
for the Mach 5 and Mach 4 nozzle configurations, 
respectively, for both the LOX and no-LOX cases. 
The corresponding value of ratio of specific heats 7 
is 1.39. 

Concluding Remarks 

A computer program, HTT, has been developed 
which calculates the thermodynamic, transport, and 
flow properties of equilibrium chemically reacting 
oxygen- enriched methane-air combustion products. 
This program is tailored to compute the various 
test stream flow properties and the mass flow re- 
quirements of the Langley 8-Foot High-Temperature 
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Tunnel for the nominal Mach numbers of 7, 5, and 4 
with and without oxygen enrichment. However, this 
computer program can be applied with only minor 
program modifications to other facilities that use 
methane-air-oxygen combustion products as the test 
medium. The option to compute the test section flow 
using flow survey data or an isentropic expansion pro- 
cedure is available in this code. The flow properties 


and operational characteristics for the expected op- 
erating envelope of the modified tunnel are presented 
in the carpet plots to assist tunnel users in preparing 
test plans. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
May 20, 1992 


Appendix A 

Chemical Composition 

The equilibrium chemical composition of the 
oxygen-methane-air combustion products is consid- 
ered the unique function of partial pressures p r and 
temperature T . These combustion products are 
thought to be thermally perfect and composed of 
4 elements, namely C, H, O, and N, and 10 chem- 
ically reacting species numbered from 1 to 10 in the 
following order: H2O, C0 2 , CO, 0 2 , H 2 , N 2 , H, O, 
OH, and NO. The six independent chemical reactions 
are expressed as 

C0 + H 2 0 = C0 2 +H 2 (Al) 

2C0 2 = 2CO + 0 2 (A2) 

H 2 + 0 2 = 20H (A3) 

H 2 = 2H (A4) 

0 2 = 20 (A5) 

0 2 4- N 2 = 2NO (A6) 

These reactions are represented, respectively, in 
terms of equilibrium constants K p j and mole num- 
bers <7 l as 


K p , 1 — 02 ^5/^1 ^3 

(A7) 

K p,2 - 

(AS) 

$ 

CO 

II 

q 

4*. 

q 

(A9) 

Ap, 4 — er 7 / cr ^ 

(A10) 

Kpfi = &&/&A 

(All) 

Kpfi = O’iq/ ^4^6 

(A12) 


The expression for K p j is given by 

log K pJ = B 1 /T + B 2 + B :i T + B4T 2 + B5T 3 


where the constants B\ to are incorporated in the 
code; these constants were obtained from reference 7 
for the temperature ranges from 360° R to 1800°R, 
1800°R to 5400° R, and 5400°R to 10 800°R. The 
units of K p3 for reactions expressed by equa- 
tions (Al), (A3), and (A6) are nondimensional, and 
those for equations (A2), (A4), and (A5) are given in 
atmospheres. 

The four elemental balance equations are 


< 7 h = 2cri + 2 o’ 5 + 0-7 + 0-9 (A 13 ) 

o*o = cf"! + 2o- 2 +0-3-1- 20-4 + o-g + 0-9 + a \ 0 (A 14 ) 

a N = 2 o- 6 + o-io (A 15 ) 

ctq = ^2 + (A16) 


The composition is calculated by solving the six in- 
dependent chemical equilibrium equations (eqs. (A7) 
to (A12)), which simultaneously represent the chem- 
ical reactions (eqs, (Al) to (A6)) and four elemental 
balance equations (eqs. (A13) to (A16)). 

For temperatures below T r = 536.67°R, only the 
major species of 04, cr 2 , 0-4, and 0's, which corre- 
spond to H 2 0, C0 2 , 0 2 , and N 2 , respectively, are 
assumed to be present in the mixture. The remain- 
ing minor species are set to zero in equations (A 13) 
to (A 16) to give 

<t h = 2 <ji (A 17) 

a O — &i + 2 <j 2 + 204 (A 18 ) 

0“]^ = 2 og (A19) 

<tq = <?2 (A20) 

The major species are directly calculated from these 
equations once the elemental composition is known. 
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Appendix B 

Thermodynamic Properties and 
Temperature 


where the values of (the heat of formation of the 
pure species) are referenced to T r — 536.67°R and 
are obtained from reference 9 and incorporated into 
the code. The entropy of the species S is given by 


Thermodynamic Properties 

The thermodynamic properties, namely specific 
heat c p at constant pressure, enthalpy h, and en- 
tropy s of the mixture at a given temperature and 
pressure, are calculated. 

The molar heat capacities of the species C p at 
constant pressure are given in terms of temperature 
as 



where Si r is the reference entropy of the species at 
T r = 536.67°R and p r = 1 atmosphere, which are 
obtained from reference 9 and are incorporated into 
the code. The value of h is calculated from 


C p<i = Ai + A 2 T + A 3 T 2 + A 4 T 3 + A b T A (Bl) 

The Cp values for the species are obtained from 
reference 9 for four different temperature ranges 
from 180°R to 720°R, 720°R to 1800°R, 1800°R 
to 5400°R, and 5400°R to 10 800°R. These values 
have been fitted with a fourth-degree polynomial to 
give the coefficients A\ to of equation (Bl). 

The value of c p is calculated using the equation 
c p = VC P t*i (B2) 


s = HtTiSi (B6) 

Temperature of Combustion Products 

The temperature of the combustion products is 
defined at a given pressure for a given fuel and ox- 
idizer and at an initial temperature of the precom- 
bustion species. An initially assumed temperature is 
iterated until convergence is obtained by using the 
following equation: 


The enthalpy of the species minus the refer- 
ence enthalpy Hi r at T r = 536.67°R, is 

r T 

Hi - H i>r = C Pji dT (B3) 

The enthalpy of the gaseous mixture h is computed 
using the equation 

h = [{Hi - Hi,) + QiA ( B4 ) 


[{Hi - H^) + QiA 
= J* j [{H j -H j r) + Qjr] (B7) 

The right-hand side refers to the precombustion 
species and the left-hand side refers to the combus- 
tion gas mixture. The initial temperature of the pre- 
combustion species is chosen to be 536.67°R. 
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Appendix C 

Transport Properties 

The viscosities of all of the species are calculated 
using the following equation (ref. 10): 


Hi = 1/1.4882 


(26.69) y^MW/T) 


(C ft) (Cl) 


is calculated using the following equation (refs. 10 
and 11): 

M ~ !&% (^2) 

where 

[i + {m/ /ij)C 2 ( M Wj/MWj) 1 / 4 ! 2 

<t>ij = -C-p (C3) 

? |2>/2[1 + (MWi/MWj)] 1 / 2 ! 

The thermal conductivity K is given by the 
Eucken equation (refs. 10 and 11) as 


The numerical value of 1 /l. 4882 is the conversion fac- 
tor from kg/m-sec to Ibm-ft-sec. Values of the col- 
lision integral Q are obtained from reference 10, in 
which fi is tabulated as a function of kT/X. Values 
of £ and A jk for the species are also obtained from 
reference 10. The viscosity of the gaseous mixture /i 


K=[c p + (5/4 )Ra]fi (C4) 

The Prandtl number Np r is calculated from 

N Pv = tic p /K (C5) 
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Table I. Typical Test Stream Gas Composition in Terms of Mole Fraction 


Conditions 

Mole fraction of species 

h 2 o 

C0 2 

02 

n 2 

Mach 7 nozzle configuration; TTC = 3560°R 

No LOX 

0.15200 

0.07599 

0.04206 

0.73000 

LOX 

.15420 

.07710 

.21000 

.55870 

Mach 5 nozzle configuration; TTM = 2350°R 

No LOX 

0.08314 

0.04157 

0.11810 

0.75720 

LOX 

.08358 

.04179 

.21000 

.66460 

Mach 4 nozzle configuration; TTM = 1640°R 

No LOX 

0.04778 

0.02389 

0.15720 

0.77110 

LOX 

.04792 

.02396 

.21000 

.71810 
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Figure 1. Operational envelopes for propulsion facilities. 
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Figure 2. Schematic of Langley 8-Foot High-Temperature Tunnel. 
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Figure 3. Triple exposure of model during its entry into test section with flow survey apparatus in stowed 
position. 
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(a) Mach 7 configuration. 



(b) Mach 4 and Mach 5 configurations. 

Figure 6. Modified nozzles for alternate Mach number capability. All linear dimensions are given in feet. 
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Figure 7. Computational sequence of HTT code with oxygen enrichment 

























Figure 7. Concluded. 
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Figure 8. Locations of flow variables for flow computations. 


















(c) Unit Reynolds and Mach numbers. 
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(d) Dynamic pressure and stagnation point heating. 
Figure 9. Concluded. 
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(a) Pressure and temperature. 
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(b) Density and velocity. 

Figure 10. Test flow characteristics for Mach 5 configuration with and without oxygen enrichment for 
TTC = 3560°R. 
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No LOX 



(c) Unit Reynolds number and mixer pressure. 



(d) Dynamic pressure and stagnation point heating. 
Figure 10. Concluded. 









(b) Density and velocity. 

Figure 11. Test flow characteristics for Mach 4 configuration with and without oxygen enrichment for 
TTC = 3560° R. 
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(c) Unit Reynolds number and mixer pressure. 
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(d) Dynamic pressure and stagnation point heating. 
Figure 11. Concluded. 
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(a) Pressure and temperature 
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(b) Density and velocity. 

Figure 12. Test stream properties for Mach 5 configuration with and without oxygen enrichment 









1000 r- 



(c) Unit Reynolds number and mixer pressure. 



(d) Dynamic pressure and stagnation point heating. 
Figure 12. Concluded. 
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(c) Unit Reynolds number and mixer pressure. 



(d) Dynamic pressure and stagnation point heating. 
Figure 13. Concluded. 
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(a) Combustor conditions. 



(b) Mixer conditions. 

Figure 15. Correlation of combustor and mixer conditions for Mach 4 operation. 
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